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MATHEMATICAL SIMULATION OF POWER CONDITIONING SYSTEMS. VOL.l; 
FINAL EEPOBT SIMUIATION OF EIME»TARY REPORT OH 

SIMOLATIOH; METHODOLOGY 


R, Prajous, J, Mazankine and J.C, Ippolit o 
Centre National de la Recherche Scientifique, Toulouse (France) 
Lab. d* Automat Ique et d' Analyse des Syst femes 


Part 1; Buck. Boost. Buck»-Boost 

1. Introduction /I* 

The purpose of this first part is to describe from a 
theoretical and conceptual standpoint the methods and algorithms 
used to simulate the elementary units buck, boost and buck-boost. 

To begin with we will discuss important general points: 

- the definition of similar converters, a definition which 
determines the generality of the results obtained and the 
effectiveness of the simulation (there is a safeguard 
against simulating the same system several times); 

- the definition of reduced parameters which characterize 
the converters and derive directly from the idea of 
similar converters; 

- the transition from reduced parameters to machine magnitudes; 

- the calculation of the corrector gain giving at the minimum 
the required "except everything" regulation of ±2Jf; 

- the general philosophy of simulations and "optimizations" 
to be carried out. 

Then we will discuss in the following order the different 
parts of the simulations to be done: 

- study of local stability; 

- investigation of a good corrector network; 

- measurement of the output impedance, input impedance 
and input-output impedance; 

- study of overall stability. 

Finally, we will end up by discussing the utilization of the 
results obtained by the simulation and we will give an example. 

2. Electric Systems and Similar Converters /2 

Let Fig. 2.1 represent 2 electric circuits and on 
which we impose the following restriction: and have the 

same number of elements of type R, L or C. At any moment t. 


* Numbers in the margin indicate pagination in the foreign text. 
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such a circuit is thus completely determined if we know the m 
currents passing through it and the m voltages at the terminals 
of the m elements. We label these currents and and 
these voltages and for and ^ respectively, with the 
subscript i varying from 1 to m. 
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^ tZ] 






\ £1 




Figure 2.1 


Definition : 

Circuit £2 is similar to circuit if and only if: 

V \ that: Vi€.Qi,...»3 

Vi e 

It is obvious that if we know the behavior of the circuit 
during the time interval [t^jt^] then we know the behavior of 
the circuit during the time Interval [k^t^^, kjptg]. 

Among other things; 

- any non-periodic steady state of circuit corresponds to 
a non-periodic steady state of circuit Cgj 

- any steady state of period T for circuit corresponds 
to a steady state of period k^ T for circuit Cgj 

- if such a state is stable for it is also stable for Cg. 



51 f" 


It Is clear that the above definition defines an equivalent 
^ ■ { ^* %»•••• J Which Includes only circuits which are 
similar to one another. 

Thus In order to analyse all the electric circuits . It Is 
. sufficient to analyse all the equivalence classes, l«e. In 
practice to analyse one clruclt per class * 

With respect to the studies to be done within the framework 
of this contract, such an idea is very important because it 
guarantees that a minimum quantity of simulations will be done 
but that nevertheless all the desired information will be present . 

It is very easy to see that two circuits and 
only similar if; 

- at two similar instants t and T they have the same 
electric structure, 

- any element R2, Lg, Cg of the circuit Cg is related to the 
corresponding element R2, L2, C2 of circuit by one of 
the following relationships: 

>S*S7’s‘i ' • 

3. Reduced Parameters Characterizing the Converters /4 

We have Just shown the Importance of the idea of a converter 
class. 

It is now fitting to numerically characterize a given class 
with the help of parameters. These parameters are necessarily 
dimensionless . 

Moreover, we have tried to give these parameters the greatest 
possible practical significance. Thus, electronic engineers 
should find them easy to use. 



Converters 

The converters we are considering may be represented as 


shown in Pig, 3.1. 


o 

Pig. 3.1 

The converter Is powered by a voltage source V^, 

The output network always consists of the output capacitor 
C (which is given an internal resistance r) and the charge, 
represented by a resistor Rg, at the terminals of which appears 
the output voltage Vg. 

Besides the above elements, the power stage Includes the ^ 
basic element consisting of the storage Inductor L through 
which passes a current 

The three converters under consideration differ from one 
another only in the way in which Vg, L and the switches are 
connected to one another . 

The circuit shown in Fig. 3.1 suffices to define the main 
magnitudes V and 1 of the power stage as well as the operating 
frequency of the converter: 




<1 



(3.1) 



Note; transformer buck-boost . For this converter it is sufficient 
to express all of the magnitudes V, I of the power stage by 
basing them in the traditional way either on the primary or 
secondary winding of the transformer. Everything we say will 
then be equally valid in this case. 


In addition, the parameters will be defined as follows; 

- starting with the steady state of the system; 

- idealizing this steady state at the output; 

- C“* , I, . C“* (the residual ripple of the output is 
ignored) . 


3*2 Parameter M, Ratio of Input and Output Voltages 
We set up the following equation; 



(3.2) 


This parameter , the physical meaning of which is very clear, 
characterizes the operating point of the system . 


Obviously we have the following; 

- KXX I o< N < 1 

- noR I 1 < N < 

- MCX-B006T I 0 < II A ♦ 

3.3 Pareuneter A, Current Ripple in the Inductor /6 

The current I^ always has one of the two forms shown 
in Pig. 3.2. 



2.T 0| T iT 

Cb) 

Fig. 3.2 
If we assume (Fig. 3.2): 

( 3 . 3 ) 

We then set up the following: 

( 3.^0 

This parameter characterizes above all the inductance 
value L: the greater L is. the smaller A is and vice versa. The 
value of A likewise enables us to know immediately If conduction 
is continuous or discontinuous. 


We simply show that: 
( 3 . 5 ) 












This parameter establishes the Importance of the "filter 
made up of the Inductor L and the output capacitor C by comparing 
Its Imaginary tuning frequency with the operating frequency. 

. The smaller » Is, the smaller the residual ripple (effect 
of r not Included) because It varies as ♦ . 

. The smaller ♦ Is, the smaller the high frequency gain of 
the power stage, effect of r not Included; It also varies 


*2 

as V . 

- *?.lnce the value of L is fixed by A, 


Infant Gharacteriaes 


for a given A, the value of output capacitor__C . 



r is the ratio between the operating frequency and the 
frequency at which the output capacitor becomes purely resistive. 

- If r is small, the quality of this capacitor is good. 

^ If r is large, the quality of this capacitor is bad. 

- In addition, T enables us to get an idea of the effect of 

r on the dynamic behavior of the power stage, the threshold 

being situated very roughly at f - 1. 

- The same Is true for the effect of r on the residual ripple. 

These last two points are summarised in Table ( 3 . 8 ). 
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Table (3.8) 


Effect of r 

r<i 

r>i 

On the dynamic behavior 
of the power stage 

negligible 

not negllbible 

On the residual ripple 

negligible 

Ripple about r 
times greate' than 



if r = 0 


^.6 Other Reduced Parameters 

In addition to the M above parameters, we can define 
other useful, but less important, parameters. 

^.6.1 Parair -^-^- -t , of the Output Current 

We set up the following: 


(3.9) 


This parameter is defined as follows: 

- for a given converter, l.e. V^, Vg, l, o, e, 

- when its output current is varied, l.e. the charging 

resistance Rg. 

The parameter J cannot be used to characterize a class of 
converters since, by definition, Ig eannot assume several values 

for the same system. 

By contrast, J Is useful In studying a given converter, as 
we shall see later on. 



3.6.2 Parameter g. Saturation of the Magnetic /9 

Cores 

For a magnetic circuit, in principle a transformer, 
we set up the following equation: 



sAturatlon 
B 


(3.10) 


We will then always take 6 = 1.33* 

3.6.3 Maximum current in the Switch 

For a real system the current in the switch 
transistor for transistors is limited to a given value. 

Without having to define a new symbol, we agree that: 


X, limitation 

L 

X. max nominal 

la 


1,33 


3.6.4 


Reduce Time Constants 


€ 


( 3 . 11 ) 


It is obvious that any time constant, in order to 
bt- reduced, must be reduced to the period T. 


We therefore set up: 



_r 

T /all 


( 3 . 12 ) 


Prom this it follows that P and the cut-off frequency 
corresponding to t are in the ratio : 


1 p 

"27ft • 


(3.13) 
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3»7 Remark It Modulator Gain 


/lO 


The definition of the similar circuits given in section 
2 can of course be applied to a converter in its entirety. 

Nevertheless, we will consider the basic circuit shown in 
Pig. 3.3 



Pig. 3.3 


Key: A) modulator D) AC return 

B) switching instants E) main return 

C) power stage 


This circuit reveals three distince parts of the converter: 

~ the power stage, 

- the modulator, 

- the correcting networks G^(p) and 02 (p). 

The similarity rules should be strictly applied to the power 
stage. 


By contrast, we can see that although the modulator also 
appears in the form of an electric circuit, it is not necessary 
that this circuit obey the similarity rule in its totality, but 
only at its output and input . 



In other words, in comparing a converter 2 to a converter 1, /II 
it suffices that; 

^ . Nr (3.111) 

In particular, the characteristic equation of the coincidence 
modulator; 


for t-t (3.15) 

e “ 

may be written aa follows; 

k v^(t) - k ^ttj - 0 for ^ 

From this it follows that the gain of the coincidence 
modulator is unchanged if we multiply the amplitude of the 
sawtooth by k and the control signal by the same factor k. 

For the sampling-delay madulator with a gain of ■ T/A, 

M 

we end up with the same conclusion. 

Aa for the magnetic modulator, it was directly defined in 
a reduced way in the "Implementation Report." 

3.8 Remark 2; Correcting Networks 

The correcting networks are defined by their transfer 
function and not by a circuit . 

For two similar converters, it th ;s si;f flees that the 
reduced time constants are identical :*nd that the gains; 

- are the same for Qj^(p) 

- are in the ratio ky/kj for Q2 (p)« 
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i*_Tranaltion from Reduced Parameters to True Parameters 
gy to Maohtne Param»ti»ya 


The probaera which we pose jn this section Is to go from 
reduced parameters to finding a class of systems to the parameters 
of a particular converter of the class. Naturally this converter 
may be the one simulated and the equations which we are going 
to give are thus also those used for implementing the simulation. 

To make this transition from reduced parameters to true 
parameters it is of course necessary to start with numerical 
values which characterize the converter in question within the 

class. To do this, it is clear that the values for the following 
must be known at the outset; 

- a voltage, 

- a current, 

- a time. 

It is equivalent and more convenient if we know the following; 

- a voltage, 

- a resistance, 

- a time. 

In, fact, the followin g equations assume that R and the 
operating period T are known . 

^.1 Equations Specific to the BUCK 
In continuous conduction (A<2): 



In discontinuous conduction (A>2); 







4.2 Equations Specific to the BOOST 
In continuous conduction (A<2 M) : 






(4.3) 


In discontinuous conduction (A>2 M) ; 



4.3 Equations Specific to the BUCK~BOOST 
In continuous conduction (A<2 (M+1)): 



(4.4) 


(4.5) 


In discontinuous conduction (A>2 (M+1)); 



(4.6) 


4.4 Equations Common to BUCK, BOOST, BUCK-BOOST 
It follows directly from the definitions given in 
section 3 that: 




P 

r ■ " " ■ ■ " 

2 7T.c.r 


(4.7) 

(4.8) 


(4.9) 
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Gain of the Correcting Amplifier ensuring the ±2$ 

Regulation Required 

A given converter Is essentially subjected to two pertur- 
bations: 

- the variations In Input voltage 

- the variations In output current I„. 

The problem here Is to calculate how these perturbations 
cause the output voltage Vg to vary as a result of their 
simultaneous action and, In particular, to calculate the minimum 
ga^ Gq of the network G^(p) ensuring In these conditions the 
± 2 % regulation required by the specifications (cf. the reports 
of the ESTEC-LAAS). 


It is useless to complicated calculations to make this 
determination. It suffices to calculate the maximum variation 
In the control voltage of the modulator (V -V ) when 
Va and Ig vary In the specified limits, then to write that: 


V - V 
e max c a>ln 


- 4 % V_ 


( 5 . 1 ) 


The problem thus becomes to determine V and V 

- . c max c min 

for each converter, and for these values it Is clear that we 

have: 


For the coincidence modulators 
and sampling-delay modulators 


‘oh 

MX 


A 

T 


V ■ 

o min ON 
min 


A 

T 


»c - - *0M If > 

»<*«<” 


( 5 . 2 ) 


For the magnetic moculator 
(cf. the Implementation 
report ) 


( 5 . 3 ) 


All of the converters to be considered operate at J « 1 in /15 
continuous conduction (of* Report of the ESTEC-LAAS Meeting of 
June 25, 1975 ^ p. 5 ) then end up by changing Into discontinuous 
conduction when J decreases. We know that: 

In discontinuous conduction, t^^ decreases when J decreases, 

- in continuous conduction, t^^ Is independent of J, 

It is then very clear that: 

** tgg occurs in continuous conduction for a value of M to 
be determined, 

- occurs in discontinuous conduction for: 
min 

* Jmln * -05. 

In these conditions very little error is committed by 
assuming the following in all cases: 


min 


5,1 Buck 

In continuous conduction we have: 



Hence t^j^ - - 0.9 T (of. section 6.1) 

max 

IMI " f ^ 


(5A) 


15 


and 


V 


«o 



From equation (5.1) we then obtain: 


®0 4%V, 




For the coincidence and 
sampling-delay modulators 




BUCK 



(5.5) 


Thus in the case of the 
simulation : 




°0 BOCK - 


(5.6) 


For the magnetic modulator equations (5.3) gives us: 


'^c ux " ^ 2T * “ 

■in 

Bln " 2T * 


Vfe see that with this modulator we cannot exceed M = 0.666. 
Under these conditions: 


V . = 0 

c min 


Hence : 


^0 4% V_ 


^0 BUCK 


- 25 .T* 


Thus, in the case of the 
simulation : 


0 .25 ^ , since it is assumed V^»Vp 

0 


9o BOCK " ** 


(5.7) 

For the magnetic modulator 

(5.8) 


/ 
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5.2 Boost 

In continuous conduction we have: 


T V -V 
OH _ 8 h 


Vrr 


«>i 

T Vg W 


H -1 - 

Hence ^)W * ^ ”5“ "I section 6.1) 


and 


T h h 

T f ■ 2 


V 4-0 
c ■in 


Equation (5.1) gives us ®o " 4%^ 


m 


®0 BOOST 



(5.9) 


Thus, in the case of the 
simulation: 


.. , o,t 

°o " ^ 


°0 BOOST " 


For the coincidence 
and sampling-delay 
modulators 


/ 

( 5 . 10 ) 


For the magnetic modulator it happens that : 


0 Bln 


T 3 I 

• ''i ■ I 5 t * " r 
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Hence : 



(S.ll) 


- 19 


“O BOOST ‘T \ 


(5.11) 

For the magnetic modulator 


Thus, for the simulation 



5.3 Buck-Boost 

In continuous conduction, we have: 


Hence 


^OM BM 


T 



♦ I 


T " IVgl ♦ |V^| M ♦ I 


( 5 . 12 ) 


Since we consider that M » 1 (cf. the Report of the EXTEC- 
LAAS meeting of June, 25, 1975, p.5): 


Bax 


T 


1 T 

I" 2 


In these conditions we obtain the same result as for the 
boost and it Is therefore not necessary to re-do the calculations 



For the simulation 


(5.13) 

For the coincidence and 

sampling-delay modulators 

(5.1^) 


BUCK-BOOST 


3*2 







For the three types of converters 



( 6 . 3 ) 

iS.H) 

( 6 . 5 ) 

( 6 . 6 ) 

( 6 . 7 ) 

( 6 . 8 ) 




The extreme values of the given A’s are compatible with the /21 
extreme values of the given M*s and with the converters; 

- operation In continuous conduction from J ■ 1 to 0.1 

- operation In discontinuous conduction from J » 0.75 to 
0.05 

(cf. ESTEC-LAAS Meeting of June 25, 1975). 

For this compatibility we used the variation equations of 
A with M and J which we will give In the section which makes 
use of the results. 


The variation ranges for the reduced parameters defined 
above are correct for each parameter taken separately . It may 
happen that some combinations are not realistic. In particular 
If A Is small (large L) and ^ is large (C very very small), It 
is possible, when F / 0, to end up with a very large value of 


r, larger than the charging resistance R_! We have therefore 

5 ' ' ' ' ~ 


eliminated any combination 


resultlnj 


resistance 


treater than R^/IO. 


6.2 Quantity of Simulations to be Done 

The preceding section shows that we have: 

- 3 values of ^ 

- ^ values of r 

- 6 values of A 

- 5 values of M for the buck and boost and 1 for the buck- 
boost . 


We can thus study: 

- 360 classes of the power stage for the buck 

- 360 classes of the power stage for the boost 

- 72 classes of the power stage for the buck-boost, l.e. 
l.e. a total of 792 power sta-^'^ classes. 
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In addition, we must consider the following: 

- 3 types of modulators, 

- 5 types of correcting networks. 

Thus, in effect, we are considering 11,880 classes of converters . 

For each class we are Interested in the following: /22 

- the local stability, 

- the overall stability, 

- the response at a step of output current, 

- the output, input and input-output Impedances. 

Thus we end up with 71,280 elementary results to be obtained . 


We give these figures to emphasise that : 

- it would be absurd to want to obtain and formulate such 
a mass of results, 

- at the limit, it would be nearly impossible to use this 
mass of results, 

- the simulation is a powerful tool of analysis: all of the 
desired results may be obtained for one well determined 
system, 

- the difficulty which exists, correlatlvely , in obtaining 
from them synthetic results. 

The rest of this section is going to explain how we avoided 
this obstacle, i.e. how we reduced the quantity of elementary 
results in order to bring it to a realistic, though high, value 
while keeping a significant value for the results . 

^ 3 Simulations and ”Optimizat ions" to be Carried Out 

With respect to the different elementary results to be 
obtained, the following general remarks can be made: 
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a) The local stability: 

- perhaps obtained using a mean calculation time; 

- the search for a local stability limit only makes 
sense for a non "optimized" system, l.e. a converter 
whose ampllfler-correcter is a simple flat gain. In 
effect, we suppose, for example, that we have determined 
that such a converter operates "best" with a correcting 
amplifier of the type G I+t^/P and with G = 5 and 

T = 0.1. It is clear that at this moment It is 
meaningless to vary G In order to determine a local 
stability limit. 

b) The overall stability: demands too much calculation / 

time . 

c) The response at a step of output current: 

- demands rather little calculation time; 

- is meaningless for a system functioning at its "best." 

d) The determination of input, output and input -output 
impedances : 

- demands enormous amount of calculation time; 

- is meaningless for a system operating at its "best." 

In addition, the reduced parameters A,t,r and M have 
different implications; 

a) r and t are genuine construction parameters . Their 
value may be chosen intentionally by the designer: 

A enable.** him to decide In what current range the con- 
verter in question is going to operate in continuous or 
discontinuous conduction or to choose the Inductance I ; 

^ then enables him to choose the size of the output 
capacitor. 

b) r relates to an imperfection, that of the output apacltor. 
The designer does not choose a value of T. Rather he 
tries to use a component of which T is as small as possible 
or as a tolerable value. 






c) M relates to an operating point. In a given application, 

M varies because varies . Thus , M has a nominal value 
and may vary by a certain percentage more or less around 
this nominal value. 

Taking Into account all of the facts which we have Just 
enumerated, we adopted the following procedure which may be 
broadly outlined as follows: 

- exhaustive study of the local stability with coincidence 
modulator, 

- simplified study of the local stability for the two other 
modulators , 

- investigation of the better correction system for all of 
the doublets (A,4>) with the conlncldence modulator and 
study of the variations caused by the extreme values of 
M and by f = 3 or 7, 

- readjustment of the correction system for the two other /2^ 
modulators and for one value of <I> and two values of A, 

- plot of Impedances and study the overall stability 
for only the optimized systems. 

This procedure is shown in more detail in the synoptic 
diagram. Table 6.1, for one unit (buck or boost type). For the 
buck-boost unit the diagram has to be modified very slightly to 
allow for the fact that only one value of M is to be considered. 
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\ <P - 0,033 

^-♦continuous conduction 
discontinuous conduction 

4 ~ 

• 1 central value of m 

P-O 

* 

Adjustment of the 
correction system 

I 

/ \ 

With sampling-delay With 
modulator magnetic 

modulator 


Table 6.1. Diagram of the simulation procedure for one type of unit. 
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7. study of Local Stablllt; 


When we consider a linear system subjected to constant 
inputs, this system possesses an equilibrium state Z® which is 
such that (using variable of state notation): 


*(t ) - 


z(t) - z® . ^ 


If the system is described by the traditional equation: 

Z . A 2 4 B U 

(with, in our case u » u® = c.^) 
it is easily shown that: 


If we set : 


B U® 


^ - z - r 


it is also easily shewn that 


Az(t) - ^(t-t ). ^(t ) , $ = continuous state (7.1) 

— o o transition matrix 


The system is asymptotically stable • when 

»(x) has all its eigenvalues less than 1 ( abs olute value) fo r 

Vx > 0. 


If ye are i 
takes the form: 


dealing with a discrete Z system, Eq. (7.1) » 


i^Z. - P" ’ 


(P = j-^'i-crete state (7.2) 
transition matrix 


We arrive at the asymptotic stability f r om the eigenval^ £ 
of 4> in the same way as in the continuous case^. 
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The llnearltj^ of the system implies that: 

- is unique , 

- the definition of the stability is valid whatever the 
value of AZq . 

If we now consider a discrete non-linear system, this 
system may be described by the autonomous non-linear recurrence 




f(0 


(7.3) 


and the equilibrium point is defined by. 

z®-£(z®), this is first order cycle (7.4) 
” of the recurrence (7.3). 

Taking the definition: 

( 7 . 5 ) 


^t is clear that: 

If A 2 ® 

■ ■ n 

Kith 

-n ~ 


(7.6) 


Eqs. (7.6) define the linear appr oximation of the system 
around Z® . 


The discrete non-linear system (7.3) is said to bo locall j^ 
stable if its linear approximation around Z ® is asymptotical ly^ 

stable . 


In contrast to linear systems: 


- Is not necessarily unique, 

- there may exist other equilibrium points than the first 
order cycles, cf. 

- ^ is infinitely small. 

The converters which we are Interested in here are discrete 
non-linear systems [ 5 ] for which we know how to determine the 
local stability [ 5 ]. 

There exists a basic difficulty in analyzing the local /28 

stability using a simulation of electric behavior. This 
difficulty is connected to the fact that the local stability is 
a purely mathematical idea in particular calling for a transition 
to the limit . 

When a linear system is unstable and when it is "initialized'' 
at a value Zo / , the absolute value of ^ Increases 

indefinitely with n. 

When we are dealing with a locally unstable non-linear 
system, the amplitude of ^ remains limited . The system 
oscillates, either at a subharmonic frequency of F = 1/T (m order 
cycles, cf. [2]), or in a quasi-stochastic manner (cf. [3-^]). 

We eliminate the non-reallstlc case where other stable first 
order cycles would exist apart from Z° (these cycles would 
correspond to non-osclllatory operating points different from the 
nominal point and which are eliminated by the electronic engineer 
for any real system by adequate safeguards). Thus, in all cases, 
the instability is characterized by the disappearance, in the 
steady state, of the strictly periodic operation of period T . 

By way of example, we can look at such oscillations in the similar 
graphs in Figs. 7.1 and 7.2 where the comparison saw-tooth <A(t) 
is given to serve as a time base. 


In principle, the oscillation is present for various 

magnitudes existing in the system. In particular, it is 

impossible for the system to be unstable and for the successive 

t-... to be eaual to one another. As a result, the oscillation 
ON 

can be observed for the command voltage Vg . 

In addition, it may also be noted that Ij^ is a more 
"sensitive” quantity (i.e. "whose variations are greater") than 
since, in a first approximation, Ij^ varies as the integral of 
Vp. In fact, in the presence cf a very low frequency oscillation, 
the relative value of Ij^ varies much more than that cf V^. 

Since the information is discrete, it is necessary to have 
a value of Ij^ and a value of per period. Since ^(t) at the 
end of a t^j^ is a discrete vector of state [1], it is clear that 
we should take I^ and at the end of the t g^^ periods , which we 
will note 

Thus the method of analyzing the local stability is based 
on the following procedure: 

a) The gain G of the error amplifier is set at the value 

G calculated in section 5. 

0 

b) The vector of state Z© of the system is set (close) to 
the equilibrium point Z° . Under the circumstances, the 
voltage at the terminals of the ideal part of the output 
capacitor is set as close as possible to the deduced 
value of Z® and 1^ is set at 10? above the deduced value 
of Z® . This 10? difference is arbitrary and simply 
results in the fact that systems which are mathematically 
stable but whose stability range around Z® is very small 
are considered to be locally unstable. 

c) The system is allowed to develop long enough so that the 
transletn behavior is eliminated: we choose to take 110 T. 

d) We take the average of and obtained for 

t going from 110 T to 150 T, i.e. ave. and ave. 
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e) We determine the maximum difference of 

Vc(toj^) with respect to ave. and ave . , i.e. 

max and AV^ max for t varying from 1^5 T to 150 T. 

f) We decide that the system is locally unstable if: 



g) If the system is stable, we Increase the gain of G to 
1.5 G and we return to (b). 

h) If it is unstable, we decrease the gain of G to 0.9 0 
until the system becomes stable again and we obtain the 
value of G denoted Gp, the limit gain of local stability . 

Remark 1 : Systems which are always stable . /31 

Some systems may remain stable even if G » or may only 
become Instable for very high gains. 

In this case, the following procedure is followed. 

An increase in G causes an increase in the residual ripple 
in Vp. We measure the gain Gj ^ where V fi (t) reaches a value of 
+A or a value of 0 and we stop the Increase in gain ac the 
following value of Gj^^: 

Although the coefficient 5 Is arbitrary, Gj^^ may be considered 
in practice as an upper limit of the gain. 



Remark 2; Joint use of and 

One might think of using only doing the convergence 

test. Some systems have a non-observable Instability for Ij^CtoN) • 
An example of this is given in the analog graph shown in Fig. 7. 3* 



ftVCK. IV® M-Oir • 



Fig. 7.3. Non-observable oscillation for 1^^. 

Remark 3 ; Limit systems of indistin ct stability^. 

It is possible to describe in an approximate manner the 
behavior of the system in question by a continuous linear 
model to which we will again refer in the following section. 

By means of this description it can be shown that some 
converters possess the following property: the phase shift caused 
by the power stage remains very close to -it for a wide frequency 











r^r/. 


range. From this It follows that In a closed loop, for a wide 
range of gain, these systems a re stable but have a damped 
o s cillatory behavior with a very small damping coefficient . 

An example of such a system Is shown In Figs. 7.1( and 7 . 5 . 

In this system the gain may be varied approximately In a ratio 

Of 1 to 300 While maintaining this behavior (only the resonance 
frequency changes). 

For a real converter or for the analog simulation the 
result Is that the system reacts to noises or to parasitic 
perturbations while oscillating constantly with a high amplitude 
With respect to the perturbations. 

We call these Indistinct stability limit" systems and, in /3 
agreement with the ESTEC (cf. Meeting of Dec. 12, 1975 ), we ~ 

will study their local stability. 
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8. Study of a Correction Network Giving Good Performance /36 

8.1 Correction Networks Considered 

Using the symbols used In Fig. 3.3, the different 
correction networks which we will consider are the following. In 
order of Increasing complexity; 

G,(pl - 0 I (0 

• ^ , advance-retard centered on 

the frequency f = 1/2 ttt. 

We have chosen to take a = this value constituting a good 
compromise between the obtained phase leading which equals 
^2® at the frequency f and the high-frequency gain which equals 
5 G, thus: 


G, (p) - G I 

G (p) . G ^ 

tTp , integrator with a 

0 at the frequency 

f = 1/2 ^t; 


(8.2) 


(8.3) 


Gj (p) ■ C 


rjp 


i» y/s *t*2 



« 


combination of the (8.^) 
two preceding equations 

with T.,<T. 

<• 1 


C|(p) - c 
Cjtp) ■ 

^ 1 ♦^p 


current loop 


(8.5) 


Examination of this list of networks Immediately leads to 
the following remarks: 

- since networks (8.3) and (8.^) possess a pure integrator 
they obviously result in a static error of 0 and should 
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only be optimized with respect to their transient response; 

- networks (8.1), (6.2) and (8.5) must necessarily be such m 

that a > a (minimum gain defined In section 5. We end 
up wlth”thls conclusion for the current loop (8-5) because 
since the mean value of the current Ij. Is 0. 1" continuous 
conductance, the current loop cannot Improve the static 
regulation of the main loop. 


g ? Method of Studying the "B est" Correction Network 

Among the factors which characterize the periormance 
If a corrected system, some play the role of constraint s In the 
sense that any system whi ch does not observe a .-Iven llmlt^s 
.elected as belnm unsatisfactory . These factors are the following: 

- the local stability; 

- the precision of regulation (cf. section 5); the 

- the Indicative response, response to a 50% step In output 
current (from Ig to Ij./2) , which must enter Into a 
determined window. This window Is shown In Fig. 8.1 
(cf. Meeting of June 25, 1975). 






The preceding constraints thus have to be completely /38 

satisfied. By contrast; 

- the output, input and input-output Impedances appear 
rather like a property of a given system for which it 
is up to the designer, the person who uses the result, 
to decide whether or not it is satisfactory; 

- the same applies for the overall stability, in certain 
cases, when the range of instability In the space of the 
initial conditions is such that the designer may decide 
that it cannot be reached or that adequate limitations 
are foreseeable. 

Thus we have adopted a method which takes into account the 
above points and which is illustrated in the flow chart shown in 
Fig. 8.2. 

The main point . with respect to this flow chart, is that 
all of the studies Involved here are not done on the hybrid 
simulation system . 

Since the simulation is only a method of analysis, such a 
task would be insurmountable. 

As a result of studies in done In another connection in 
this area of the LAAS on the theoretical level, we have ^ 
approximate continuous linear models of the three converters 
buck, boost and buck-boost in continuous conduction and 
discontinuous conduction [6]. The models, although approximate, 
give a good Indication of the behavior of the system in 
question by giving the gain and phase diagrams of the approximate 
transmlttances ; 
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I Predetermination of 
the parameters of the 
system 


Adjustment of 
parameters 


Can the network 
N*N give the 
stability, the 
precision, the 
indicative re- 
sponse? 


No satisfactory 
solution 


Result good? 


Adjustment 

possible 


Study of the overall 
stability 


Overall stability 
satisfactory 


FIk. 8.2. Procedure for studying 
the "best" correction network. 


Determination of 
impedances 








from specific values of A,$,f and M and from the machine values 
Vg = 0.^; F = 10; Rg = 1 (see, for example, the diagrams in 

Figs. 7.4 and 7.5). 

The following parts of the flow chart shown in Fig. 8.2 

were done using thse models: 

- what is the simplest network which can be adapted? 

- predetermination of the parameters of this network. 

The simulation is used to: 

- adjust the parameters of the network taking into account 
the exact behavior of the system; 

- determine the exact impedances; 

- give the overall stability, which the approximate model 
cannot give . 

Thus it was possible to carry out this basic phase of the 
study of elementary units only by the Joint use of a powerfuj^ 
but approximate theoretical method on the level of synthesis, 
and a method of exact analysis using simulation . 

9. Determination of Output, Input and Input-Out put Impedance /]n 

The notion of Impedance is a notion purely related to con- 
tinuous linear circuits. 

The converters in question are discrete non-linear systems. 

This implies the following obvious consequences: 

- the notion of Impedance in the case of a converter is 
necessarily an approximate notion; 

- the corresponding measurement can be made only at a small 

signal ; 

- the significance of the values obtained increases as the 
frequency in question decreases and their significance 
decreases as one approaches the frequency F/2 which 
represents a limit (SHANNON theorem). 
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9«1 Theoretical Method Used 


Let us assume, therefore, Fig. 9.1, an impedance dipole 
Z(p) represented symbolically as follows: 

, Al(p> 

|| 2(p) f AV(p) 

Fig. 9.1. Elementary dipole. 

The problem which we have to solve is one of determining 
experimentally (in some way) the function |Z(,Ju))| for a certain 
interval of variation of the angular frequency to. 

The difficulty basically stems from the fact that for the 
physical system (or for the analog simulation, which amounts to 
the same thing) we have temporary values Av(t) and Ai(t). 

The most obvious method consists in applying to the dipole 
a signal Av (or Al) of sinusoidal form Aslnwt and measuring 
in a steady state the absolute value of the response Al (or Av) 
and setting up the ratio of the absolute values. 

This extremely traditional method certainly leads to the 7^2 
result but in our case It has the following drawbacks: 

- since the measurement is made for a small signal, the 
noise present in the machine is not negligible and we 
cannot be content with a measurement of the peak value; 
the data obtained must be processes, for example a 
correlation process, and this treatment must be carried 
out for each value of w, which extends the simulation 
time ; 

- when (i) is small, the measurement lime required for one 
vlaue of 10 is long, approximately 2 x 2ir/aj. 


By way of example, it took us about 35 minutes to plot a 
curve |z(Jw)| using only 20 values of w. 

Therefore we used the method which we are now going to 
describe and which is based on the property that Z(Jto) is the 
FOURIER transformation of t(t) . 

Therefore let : 


E(t) t £'^ [a(p)] 


(9.1) 


We break Z(p) down into simple elements: 


Z (P) • — T~ * * ' 

P -«1 P -»2 


(9.2) 


assuming that there are no multiple poles nor zero poles and 
that the number of zeros of Z(p) is less than the number of 
poles . 

We then have: 

•|t 

ft + ... (9 3) 

The FOURIER transformation F(o)) of z(t) is: 




p(ua) s / i(t.) Q. dt 


(9.i^) 


This gives us: 


r(bl) 




i41 


dt 


^ • • • 


(9.5) 





pN) 


N 1 — «-6) ai 

-v15>L Jo Vl-L Jo 


Hence : 


Aj »2 

■ 5 * 57 * 35 : 5 ;* •■ 


(9.7) 


Comparison of Eq. (9.2) and Eq. (9.7) clearly shows that: 


” Z(j*W 


(9.8) 


.r nhe method Is thus to apply a D IRA£ 

H.ni-n function A1 = I^ 6(t), to measure the output Av 

and to perform the In tegral (9.^) off-li ng.* 

Theoretically, a slnele simulation giving z(t) Is sufficient 
to obtain the entire curve |Z(3m)|, which represents a great 
improvement with respect to the preceding method. 

9 2 . a.ntetion of the Preceding Method to_0aLCa51 

„ , Determination of the ^*e F, met Ion Respons_e 

It Is clear that It Is not possible to directly 
apply a DIRAC delta function to the dipole In question. 

one possibility Is to apply an Impulse of finite amplitude 
„ for a sufficiently short time t, Fig. 9.2, and 
impulse with a DIRAC delta function of area Ht . 

„ Lst be compatible with the possibilities of the ol-u t 1 

ouestlon. Moreover, since t Is small , the 

system Is low and the corresponding response Ht s(t) 

amplitude . 


Fig. 9.2 


Another method can be used. 

The behavior of the circuit Is determined by the following 


equation of state 


X - wc + Bt 


V - CX 


(9.9) 


If i is a DIRAC delta function i = we have; 


Thus : 


X(O') ■ 0 
x(0*) - M 


X(O^) - • ^ B1 


v(t) - C • ^ B - 


(9.10) 


The impulse respo n se being sought z(t) is thus only t ^ 

of the c irndt In a free^te but conveniently set ^ 

the value BI^. 

For a real circuit, thlo adequate Initial setting of 
can generally be determined very simply by examining the circuit 




We give to Iq the maximum value compatible for with the 
assumption of small signals. 

9.2.2 Discrete Evaluation and Integration for a 

Finite Time 

In the hybrid simulation we have infact only one 
value of V per period T of the converter, given the discrete 
nature of the behavior of these systems. 

In addition, v is only measured for a finite time kT. 

Eq. (9.4) becomes: 

r’M 

As a result, modified Eq. (9.6) gives: 


1 


S(t) • dt 


(9.11) 


F'(*« 


r*(rt) - fW 


A. r (a.-JW»kT "1 Aj r 1 

’j ■ ■*]♦••• 

aj-j*o[_* aj-j«> J 


( 9 . 12 ) 


(9.13) 


error term 


In the discrete form we in fact calculate: 


F"M - 


H . .. 

n«i n 


-jOnT 


(9.14) 


By way of testing, we apply this method to a circuit of the 
second order wired on the analog part and such that: 




"t - a;! 


.1 j 


-.1, 


liocs lisu w.-. :i:L-.iLa. 


- continuous impedance is 0.1 

- resonance of 7 dB at frequency 1. 

The result is shown in Pig. 9.3 and one can see that it is 
very satisfactory . The oscillatory appearance of the curve at 
the low frequencies is explained very well by the form of the 
error term of Eq. (9.13) which, at the given kT, is sinusoidal 
in terms of to with an amplitude which increases when to decreases. 

In practice, this error is not detrimental because it manifests 
itself in the zones where the behavior of the system is the 
best known. 

Circuit of the Second Order /^6 






.ttm Mi. 


•.4 


tJ 
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Fig. 9.3. Example of Impedance obtained by FOURIER 
transformation . 


In order to make a comparison with the first method, the 
obtaining of such a curve requires 90 seconds (essentially 90 


^5 


.'1 




Ja-"' 




seconds of processing on our CII-MITRA 15 computer) for 200 
values of u) (thus for a very detailed curve). 


Since the main portion of the time is the processing time, 
in order to Improve the results while decreasing the noise, we 
in fact make 10 succ essive firings in acquiring the z(t) values 
(which requires only 30 seconds) and we take the average of the 
samples of the same order before actually processing the data. 

9«3 Output Impedance 

Let us consider the output network of the converters 
(Fig. 9.^) to which a DIRAC impulse I^fiCt) is applied. 


Ay 



V(t) 




^^8* 9«^. Application of a DIRAC Impulse of current 
to the output network. 


Immediately we see that the effect of applying the Impulse 
at t = 0 is to create at the terminals of C a discontinuity of 
the first kind of value 


c (9.15) 

On the input side the network is either in the air or 
connected to the Inductor L. In both cases this connection 
does not occur. 



We therefore obtain s„(t) by the diagram shown In Fig. 9.5. 

B 



Fig. 9 . 5 . Obtaining In practice the output Impulse 
response . 


The result Is given In the reduced form; 


1*. (iwol 

R_ 


(9.16) 


In a first approximation, the output voltage may be 
broken down Into the sum of a continuous value and a residual 
ripple whose form does not vary at small signals. It Is for 
this reason that It is sufficient to measure AV^ per period to 
obtain directly the variation in the mean value of Vg. 


9.^ Input-Output Transmittance 

By definition, this Is the following quantity: 




(p) 


(9.17) 


This Is a dimensionless quantity for which we must speak of 
transmittance and not Impedance. 



Since the measurement of this quantity requires the 

application of a DIRAC Impulse on the output, is determined 

at the same time as Z_. It suffices to measure AI. and AV„ 

■ ” — s A o 

after applying the step V^. 

The determination of the mean AI^^ is not as obvious as 
that of the mean AVg. Fig. 9.6 shows what I^ can look like 
during the course of a steady state period and a transient 
state period. 



Fig. 9.6. Variations in I^. 

We see that it is not sufficient to know AI„ in order to 

M 

know AI^ ave . 


The following procedure is then used: an auxiliary 
integrator reset to zero at the beginning of each period gives 
at the end of the period the quantity /'*' i.(t) dt - t i ave. 

•'O A A 

and I^ ave. becomes directly acessible on the analog wiring 
no matter what the form of 1^^, l.e. whatever the converter. 


Remark : 


in continuous conduction . 


Continuous T^g, l.e. Tgg(O) may be calculated in a very 
simple manner. 



In the steady state we have; 


V. I, Doy - I_ 

A A ^ s s ave. 


(9.18) 


way V 


assumed that the j;ain is sufficient /50 

S A ... . 


We therefore have: 


such that V s c.^ ) • 



(9.19) 


which gives us a very simple way of verifying the curve |T„e(Ju))| 
obtained . 


9.5 Input Impedance 

Given the electric nature of the system connected to 
the input terminals, in order to measure the input impedance it 
is clear that it is necessary this time to apply a voltage AV ^ 
and measure the variation AI ,^ ave resulting from this . 

We thus determine an input admittance . 

The correctness of the above reasoning is verified by the 
fact that the input admittance decreases with the frequency, 
which must necessarily be the case in order for the method 
discussed In section 9.1 to be applicable (number of zeros less 
than the number of poles). 

Depending on the moment considered within the period, the 
source is connected to the Inductor L or it is not connected 
(Pig. 9.7). 


^9 






Fig 
We define: 


.9.7. Connections of to the converter. 


time in which L is connected (g. 20 ) 


/5I 


K 


suppose that we add to a DIRAC Impulse V^6(t). This 

rrinsT’v it does not have any meaning 
procedure is purely imaginary, it ..^vprter 

as such because of the discrete nature of e 

Nevertheless, we can deduce from it an 

1-hp svstem If L is connected to apply ng 

on the system, j- ^ , j j t and there is no 

causes a discontinuity of the first kind on ^ 

"stantaneous effect on the other parts of the circuit. The 

equivalent discontinuity on '=°- 


L 


.-i K 


C9.21) 


For all of the converters the Impuse response at the input 
is obtained by means of the circuit shown in Pig. 9.8. 

Remark 1 : Values of K_ 

The coefficient K defined in Eq. C9.20) equals: 


50 




Fig. 9.8. Obtaining the input Impulse response. 


K a tQjj/T for the buck and buck-boost in the two conduction 
modes 

K = 1 for the boost in continuous conduction 

V _ t /T for the boost in discontinuous conduction. 

ON OFF 

The following values are easily calculated; 

K-H Buck in continuous conduction (9»22) 

K-i Boost in continuous conduction (9.23) 
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Buck-boost in continuous conduction (9.2^1) 

in discontinuous conduction (9.25) 
Remark 2 ; Presentation of the result 

Although we determined an input admittance Y^(Jo)), we in 
fact plot the input Impedance z^(j^ . # presented in the 

reduced form; 




Remark 3 ; Value oT In continuous conduction . 

We have : ave 

(We assume that the gain is sufficient so that and I, 
are nearly constant . ) 

Thus : 


dV. 


noy.dV^ ♦ V^.d Boy • 0 

V** _ 


(9.27) 


dl^noy 


Ml^moy 


Key: moy = ave 


Thus we have : 


z^(0) 

A 


( 9 . 28 ) 


10. Study of the Overall Stability 

For this section we will refer to references [1-2-3-4-71. 


10.1 Definition 

Let us assume a discrete non-llnear system described 
in the space of state by the recurrence: 




- f(Z ) 

— "Tl 


( 10 . 1 ) 


A cycle C of order m of the recurrence 
definition a set of m vectors: c -Jz,j, 5c2'*** 


(10.1) Is by 

,z X all consequent 

-on \ 




In other words; 


If 


( 10 . 3 ) 


C la an equilibrium set of the recurrence (10,1) . 


C is a stable cycle If the matrix : 


R ^ 

®C d8 


5«€- c 


(10. H) 


at all eigenvalues less than 1 In absolute terms. 


There exists another type of equilibrium set called 
"stochastic behavior zone” which we only mention In the frame- 
work of this report. Let D be the closed space 'f the admlssable 
Initial conditions for the system described by Eq. (10.1). 


By definition, the system described by Eg. (10.1) has 
an overall stability If the only stable equilibrium set Included 
In D Is a cycle of the first order . 


/ 5 ^ 


This cycle of the first order corresponds to the normal 
operation In a steady state of the system. Any other stable 
equilibrium set would correspond either to a subharmonic 
oscillation or to a pseudo-stochastic oscillation. We have 
already given examples of such oscillations In section 7. 


Corollary ; If the system (10.1) Is on the whole stable 
within the range D, any state set to any value whatsoever of 
Zq D converges towards a stable cycle of the first order. 


10.2 Method Used 




Very elaborate methods exist for determining what 

5 
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the equilibrium sets are for a given recurrence, stable or 
unstable, and for plotting the limits of the stability ranges 
of the stable sets (cf. references mentioned at the beginning). 

These methods, which are especially designed for making a 
detained analysis of the properties of a given recurrence, are 
not used In our case for the following reasons: 

- they are essentially algorithms which can be used In 
numerical calculation, 

- they require much calculation time, 

the precise plotting of the boundaries of the stability 
ranges Is only of little use for the user of the system. 

As a result, we used a method based on the following points; 

a) In the case of the converters In question the recurrence 
(10.1) Is at least of the second order, since the power stage 
Itself Is of the second order. Of the correction networks given 
In section 8.1, the most complex is network (8.il) which is of 
the second order. The recurrence corresponding to the system 
corrected in this way would thus be of the fourth order . 

Analysis of this system would present very great difficulties, 
both in terms of calculation time and result evaluation time. 


As a result. In agreement with the ESTEC (cf. ESTEC-LAAS 
Meeting of Jan. 12, 1976), we decided to limit ourselves to 
systems of the third order. More specifically, we will study 
the systems of the fourth order In an approximate manner, as 
If they were third order systems. This approximation Is 
perfectly Justified. In effect, the pole p = -/^/t present In 
networks (8.2) and (8.^) is situated beyond the equivalent cut- 
off frequency of the system. 


( 

! 

I 

j 

In summary : | 

- the "optimized" converters provided with networks (8.1) ' 

and (8.2) will be analyzed as second order systems; ; 

- those provided with networks (8.3) » (8.^1) and (8.5) will | 

be analyzed as third order systems, the variable of state ! 

coming in addition to those of the power stage which is the j 

output of the main integrator of the network. ' 

b) The discrete vector of state in question is at the 

beginning of t^j^ in the case of second order systems 
and in the case of third order systems. 

Li'-J 

For each variable of state x we define a maximum value and 
a minimum value and 


Vie consider that the range d of the admissable initial 
conditions is the parallel piped limited by the different maxima 
and minima . 

D is discretely evaluated by giving to each variable of 

I - “ *mln . *max " *nln, f » 

state X t he values ♦ 3 )# 5 

The domain D is thus made discrete by I 6 points in the case 
of the second order and by 6^ points in the case of the thrid 
order . 

c) Then we set the system to each of the values fe D | 

and which we have Just defined, and we observe if the system 
converges or not towards the normal steady state (first order 
cycle). To do this we use exactly the same algorithm as for the 
]ocal stability. The only slight difference is that the con- 
vergence is not tested at 5? by at lOJf or even more, because /56 
the quality of the analog wiring used for the overall stability 

I 



Is Inferior to that of the wiring used for the local stability 
(In this regard, cf, the first part of volume 2). 

10. 10,3 Presentation of Results 

Pig. 10.1 shows an example of how the results are 
presented. We use a perspective representation of the domain 
D which has been divided Into discrete units 


A 

BUCK,STABILITE GL08ALE. M»0 .67 
GAt1-e .O0DELTAS0 .33PHI«0 .10 

c=l.5 


B Plon ZsZiNin 
• » • • 

• • • • 

• • • • 


Plan ehoifti pour 
f'ordraJL 


Man 


lb 

Mill 


■ • • ■ 


- • 

- '• 

• - 


• • • 


VftMln Vonmji 

Pian ZsZmaii 




USHAX -0.S7 
USHXN -0.00 
ILHAX »O.90 
ILHIN *0.00 



Pig. 10.1 Example of how the overall stability Is 
represented. 

Key; A) overall stability c) plane chosen for the 

B) plane second order 


Each point of the domain D is represented by an alphanumeric 

character with a precise meaning; 

- the dot means that this setting has not been analyzed; 

- the minus sign means that the corresponding setting 

converged towards the stable first order cycle; 

- the zero "0" means that the corresponding setting did 
not converge towards the first order cycle (instability). 

The symbols and "0" were choses so as to be very 
easily distinguishable by sight. 


Thus we have overall stability if and only if all of domajji 
D contains only minus signs 

To be sure, strictly speaking we are not absolutely certain 
because the discrete evaluation of the domain D is not very 
detailed, but nevertheless there is very littl e chance that any 
instability present will not be detected. 

For the second order systems the results are presented in 
the same way as for the third order systems, but only one of the 
four planes corresponding to the discrete evaluation of Z is 
used, the other three then being filled with points. 

The coordinates of the parallel piped D are clearly 
indicated and for more clarity we have discussed instead of 
\ (since Vg « Vj^). 

Remark : Systems close to the local i nstability. 

- Some of these systems are locally stable and globally 
unstable. In this case, it is very common that the stability 
range of the first cycle has a very restricted surface area (or 
volume), but nevertheless its form is very complicated allowing 
it to attain very elongated zones of D of the stable first order 

cycle . 


- In addition, such systems generally meet with a high 
deeree of sensitivity and, a fortiori, a sensitivity 

— ^ which Incerases as k Increases. 


In physical terms, this means that a small variation 
is going to result, for example, in a substantial variation 

''2n+50- 

- In the case of the analog simulation, these two facts 
combine along with the imperfections and noises to lead to the 
following phenomenon: the same system, set several times in 
succession to the same value, each time follows a different 
development . 


It may happen, for example, that an initial condition 
taken out of the domain of stability of the first cycle, at the 
end of a certain number of periods may lead to a discrete 
vector of state penetrating into the domain of stability and 
remaining there. 


As a result, when we are faced with a domain D simultaneously 
conta i ning and ”0*s”, the only valid information is that 

the system is locally stable and globally unstable. It is 
impossible to get any indication of the topography of the 

and ”0*s” in D , except for certain cases which will be 
pointed out when we discuss the actual results. 


The example shown in Fig. 10.1 is precisely in this case. 

The example shown in Figs. 10.2 and 10.3 show the same system 
with a larger and smaller gain respectively of the error amplifier. 
Fig. 10.2 clearly shows that there are no longer initial 
conditions leading to stability and Fig. 10,3 shows that the 
system is globally stable. 
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USIHIN -t.97 
U8H1N «e.M 
ILNAX ••.9C 
ILHIN -9. 89 



Pigs. 10.2 and 10.3 (next page). The system shown in 
Pig. 10.1 with different gains. 


Key: A) overall stability 


59 


BUCK'.STABILITE CLOBALE. M>8 .67 
CAn»0 .OBDELTA-e . 33PH I 30 . 1 e 

I 

c-e.5 


• •* 


• • • 


t • • 


• • • 0 


0 0 0 


UtNAX -0.07 
USniN - 0.00 
tLHAK - 0.00 
ILHtN -0.00 



Fig. 10.3 


11. Utilization of the Results 

In this section we will use the following notation: 

subscript S for the variables relating to the simulation; 
subscript V (true) for the variables of the real 
converter in question. 


^ Trans ition from the Numeric Values of the 

Curves and Hard-Copies to True Valueti 

It is now assumed that we have curves or hard-copies 
corresponding to the values of A, «•, f and M of the real 
converter which we Intend to "optimize." 

The problem here is to go from the numeric values of the 
simulation to true numeric values. 


11.1.1 Gain 


all of tho Is Involved in 

results. When this amplifier does not have any 

correction function, l.e. when G Co) - G fho 4 ^ 

Tn - G, there Is no amblKultv 

h" :::.T:r: :r:: 


class‘the"irL!'"“T'“"“'“’* sam, 

Class their loop gains are Identical. This: 


°S*Ss*®PS “ V®MV*®PV 

with Gj^ = modulator gain 


( 11 . 1 ) 


<st 


'ON 


dV 


Gp = power stage gain ^ 


<Jt. 


ON 


For the coincidence modulators and sampllng-delav 
-dulators. G„ Is of the form T/A. We therefore ohtaL: 


“s 


-a !» /””s ) ■ 
*s \ “oh/, \ ; 


(11.2) 


The quantity T dV^/dt^^ Is proportional to V„. 


Thus : 


88 


8V 


* »s ®» 


(11.3) 


Hence : 


Ay V. 


88 


r, 07 i 


0,4 


SV 


Thus, in order to go from the gain read on the Gg curves to the 
ture gain, we have; 


®v 



for the coincidence and 
sampling-delay (11.5) 

modulators 


The operation of the magnetic modulator, such as it has 
been defined, leads to a constant quantity * 5 ^ T which is 
equal to 2 / 3 . in the case of the magnetic moculator is 
thus proportional to T/V^ = T/Vg since it is assumed that = 

''s- 

Eq. (11.1) then becomes: 



( 11 . 6 ) 


As above, the quantity 
Vg and we obtain: 


T dVg/dTQjj is proportional to 


G 


SS 


S* V, 


S6 



Thus : 


for the magnetic 
modulator 


(11.7) 


The same line of reasoning is carried out for the gain of 
the current loop, but Eq. (11. M becomes: 


'ACV 


Ag ^ ^ * *SV V '' 


"sv 


8S 


SV *88 


( 11 . 8 ) 
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since we have assumed R 


1, we obtain; 


SS 


for the coincidence and 
sampling-delay (11.9) 

modulators 

For the magnetic modulator, Eq. (11,6) becomes: 

( 11 . 10 ) 

The quantity T dl^/dt^j^ Is proportional to V^/Rg, hence: 

G ^ -L. . c 1 

"ss ^ ^7 (11.11) 



Thus: 


°ACV " °ACS*"SV 

for the magnetic (11.12) 

modulator 

11.1.2 Time Constants. Poles. Zeros 

This time we will consider in a correction 
network a term of the form (1 + tp) located either In the 
numerator (the zero case) or In the denominator (the pole case). 

For all of th e elements It Is clearly Indicated “pole 
(or ze ro) at F/x. ” Everything Is produced to the operating 
frequency F and. In accordance with the definition given In 
section (3.6.i<), the corresponding reduced time constant Is: 


j 


.. 


(11.13) 




1 


As a result: 


'tr' - 

V 




2n 


V 

2 11 


( 11 . 1 ^) 


11.1.3 Curves as a Func tion of xlme 
These are curves of response to a Jump in 

output current . 

They may be represented in two forms; 

a) AVc,(t). We use the variable expressed in crdinat 
values (AVg) by means of the followinr, equation. 
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1 

1 K V m 

1 -*^sv 


0,4 

b) Al^Ct). We have 





L- 


max 


cn. 15 ) 


(11.16) 


(By I is meant the value of Ij, at the end of t^^,.) 

for the seale alonr, the it ^ 1 s_^’.e^lv_div Ided 

Into the number of sampl e s, i.e. _t he number , of. rejirdil • 

11.1.^ Curves r?^ « Function of Frequency 

a) Output and Input Impedances 

since the ordinate of theae curves Is divided 

into values of IZ/Rgl. “ 


Kl- 1^1 "" 


(]1.17) 


6U 











b) Input-Output Transmittance /6A 

Direct utilization, the input-output trans- 
mittance being dimensionless. 

As for the scale along the abscissa, we have 



(11.18) 


11.1.5 Overall Stability Diagrams 

As above, we go from a quantity Xg to a 
quantity Xy by means of a proportionality relationship; 


Scale V : v_, 

s sv 



(11.19) 


Scale ij. : - ^Ls 


1,331 


LV BOX 


LS nax lu 


( 11 . 20 ) 


Scale z * 


max 

2 

S max lu 


[note: lu = read] 


( 11 . 21 ) 


11.2 Change in Operating Point of a Given Converter 

Let us assume a given converter, l.e. operating at 
a nearly constant voltage Vg as a result of closed loop 
operation and with the elements L, C and r with given values. 

The operating point of this converter changes when 
varies and/or when Rg (or Ig) varies. During such an operating 
point change, the converter changes class. Its behavior then 
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corresponds to curves other than those used previously. 


During a change In operating point, (f> and r are constant , 
but M and/or A vary . 

Below we are going to give the equations allowing us to 
calculate these variations. 


Consider a converter for which all of the values of the 


elements have been calculated for nominal values A 


nominal * 


M J 

nominal* nominal* 


If we use this converter with a different 
M (different V^) and a different J (different I 5 ), the equations 
which follow give the value of A for the new operating point . 


11.2.1 Buck (0 < M < 1) 

1. Nominal values in continuous conduction 


This is the case if 


id <.2 

nominal ^ ^ 


A . A nom 1-M 

nominal* J * 1-M 


nominal 


if 


( 11 . 22 ) 


The operating point in question is in continuous conduction. 


A 1/2 A nominal , *^nom 1-M 
“ r J ' 1-M . 


if A >2 


(11.23) 


nominal 

The operating point in question is in discontinuous conduction 


2. Nominal values in discontinuous conduction. 


This is the case if 


A , , > 2 

nominal 




A - A 


s! 


^nom 1-M If A> 2 (11. 2M) 


noalnal ” J 1-N. 


noAinal 


The operating point in question Is in discontinuous conduction. 


J 

A nominal nom 1~W 


• ■ I • 


J 1-M 


nominal 


If A ^ 2 (11.25) 


The operating point in question is in continuous conduction 


11.2.2 Boost (MM) 

1. Nominal values in continuous conduction. 
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This is the case if; 


A ^ 2 N 

nominal nominal 


A . A 


nominal 


J /'h , , 

nom I nominal | 

• J *V M / ‘ M 


M-1 if A < 2M (11.26) 


nominal 


-1 


The operating point in question is in continuous conduction. 


f 


A-y 2A 


3 M 

nom nominal M-t 


nominal* J 


H H ,-l 

nominal 


if A> 2H (11.27) 


The corresponding operating point is in discontinuous conduction. 


2. Nominal values in discontinuous conduction. 


This is the case if 


^ nominal ^ ^ **nomlnal 


07 




K A J — r 

^ “ noain«I V ^ •* *V\ottln»l 


if 


\o. op<‘ rating 


point in question ts in dlooontlnuous conduction 


fjc. ISiS^ ir A<» (11. ? 9 ) 

■ — t; .1 2 H 


nOAlMl 


operstlns Point In question Is In continuous conduction. 

11 , 2 , 1 nu ^c.Bocst (M>0.) 

1. Nominal valuoa In contlnuoua conduotlon. 

Thla la tho caso If: 1 A — • ***nc^nai*^ 




A <.2 OHl) (11 oO) 


nominal J 


The opcPStlnp point In question Is In continuous conduction. 

Tno oporstln. POlnt In question Is In dlsconllnuous conduction 
2. Nomlnsl v.slues In discontinuous conduction 
This Is the esse If ^ J 


A . A 


nominal 


If A >2 (H*n ( 11 * 3 *^ 


The operating point in question 


is in discontinuous conduction. 


» Aniln.1 !no. J_ If 

^ 5 J l+M 


The operating point In question Is In continuous conduction. 

1 1 . q Change of Freque n cy for a Olven Convert er 

The given converter Is defined as In section 11.2. 

If this converter has been calculated for the nominal values 

. * r and F It changes class If we change Its 

^nom* *nom* ^nom nom 

operating frequency expressed as F. 




We obtain: 




noa 


‘^noa 

r 


(11.34) 


r -p 


noa r. 


noa 


(11.35) 


A -A ^ 
noa r 


A mA 




if we go from continuous con- (11.36) 

duction to continuous conduction 


if we go from discontinuous con- (11.37) 
duction to discontinuous conduction. 


If we change the mode of conduction as a result of the 
change in F, we use the following procedure: 


■ nom 


nom 




Equation of 
mode 1 


Transifion between 
mode 1 and mode 2 





mode 2 

Equation of mode 2 using 
the limit values instead 
of the nomlmal values 


11.4 Method of Utilization 

In this section Ke do not claim to give the only 
method for using the simulation results. 


In effect, any designer famllar with the results will be 
liable to use them In his own Individual way. 


Our Intention here Is only to provide the following: 

- a certain number of remarks of a general character 

^ at ions • incrsaslng 
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which may guide the user in his 


safety margins, etc; 
- a utilization method 
used in the greatest 


which should be capr^-ble of being 
number of cases. 


11.4.1 General Remarks on the E f fects o X 
Variations in Parameters 


The data which follow are approxima,^ and may 
turn out to be In part Inexact for some systems. Nevertheless. 
In the majority of cases they are satlsfscolry. 


a) The converters In the discontinuous mode have their 
gain In high frequency which decreases 20 dB/decade. 

The boosts in continuous conduction have their gain In 


high frequency which decreases 20 dB/decade. 

The bucks and buck-boosts In continuous 
their gain In high frequency which decreases 


production have 
40 dB/decade. 
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b) The high-frequency gain; 

- varies as 1 /M 

2 

- varies as 

- practically does not vary with A for the continuous mode 

- varies as 1 /A for the discontinuous mode. 

As a result, If we want to keep the same pass band equivalent 
for a given system , we must : 

- use the gain 5^ 

- use the eain > / **optiD ^ ^ 

°optla \9 ) 

- use when A varies, in the continuous mode /70 

- use ^^^opt discontinuous mode 

- in the case of the current loop, vary and G as 
Indicated above. 


c) If we want to Increase the safety margin of a system by 
making Its equivalent pass band go from to 

must : 

- use the *^Qptj_jj,/x the discontinuous mode and for the 
boost in the continuous mode. 


- use the gain 

the continuous mode, 

- use G 


for the bucks and buck-boosts in 


optlm 


/X and 0^ 

OptlB ! 


d) The presence of T / 0 ; 

- Is unfavorable for the discontinuous mode and for the 
boost In continuous conduction, 

- Is favorable for the buck and buck-boosts In continuous 
conduction. 
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_ is negligible If T < equivalent pass 

h P Pronos n- « UtllUatlonJtett^ 

The saln_E0lnt racf.Tercrired 

optimisation of the system s In^gues^^ eonverter^ ^»TJ!l»Llh^ 

above, that 

nperatl np; point . 

K Investigate what the worst operating 

Therefore we shoul stability and correct the 

conditions are with conversely, «e determine 

system under these condition. . the lowest 

what are the operating cond „„tput stage and we decide 

rerr^:trrhrrrs: Eduction in performance la ac 

ceptahle . 

Taking into account what we have said In section ll.^.l. 
the high-frequency gain of the power stage 

- M is minimum and 

- & is minimum. 

relating the variations of A to thos 
The various equation increases. mhP worst case 

- ” ^%'Trt^converter delivers Its maximum current. 

The opposite case thus corresponds to the following: 

- M is maximum and 

- j is minimum. 


12 . Utilization Examples 


The following examples are examples relating to boos _t 
converters and the reader Is thus requested to refer to part 
1 of volume 3* 



Conduction 

Let us suppose a boost converter possessing the 


following characteristics; 

4 - ji ^ 40 • 

Vg - 40 V Ig - » ^aax thus "s ^ 

V. may vary from 30 V to 20 V 

F = 20 kHz thus T = 50.10 
We want the converter always to operate In discontinuous 
conduction and we expect to use It from Ig = 1 A to Ig = 
0.05 A. We can be content with a theoretical residual 
ripple on the order of 1 % or even greater, such that 
C is small. We use a coincidence modulator with A = 5 V 


The parameter M may vary from 40/30 to 40/20, thus 


1,33^ M < 2 


( 12 . 2 ) 


Eq. (11.28) easily shows that, for a boost In discontinuous 
conduction, the quantity A/M decreases when M Increases. Now 
we are In discontinuous conduction If A/M > 2. Thus as M 
increases we come closer and closer to continuous conduction. 

As a result, we begin the pr o portioning at M maximum and ,^ 



Ke Should have A > 2 i.,. , , ^ 

us we elect to operate with A = 5 when M = 2. 

I'iViCs'a; r-'®’ ^ 

-i-.jj always at Ig max: 


A - 5 1/— 2- 0.33 , 

\/i»33 • T — - 3,5 


valuel: *° Proportion the converter for the ncni.al 


*'no« - »'33 


^no» " 3.5 


^ - 0,1 (the relative residual ripple is 
very roughly equal to <p^) 


(12.3) 


|»g noo - 40_f^ 


Eq. (^.4) gives us; 


L . 2i33 2 

1.33 3,5 


1.33 1 c2 * ^0.50.10 


1* ■ 80 1^ H 


Eq. (A. 8) gives us; 


47t^.0,l^.{2.10^)^80.10‘* 


(12.4) 


c - 8O1 


(12.5) 


The capacitor we have becomes purely resistive at 02 


kHz . 



Thus we have ; 


" * (12.6) 

17i_e power stage Is at Its gain maximum for M = and 


It Is for these values that we are now going to determine 
the correction networks. 


We begin by analyzing the local stability curves for which 
we have the value A = 3.8, sufficiently close to that which we 
are Interested in. We find that for the values ^=0.1, 

1 and A - 3.8 it is Impossible to obtain the local stability 
for the gain leading to a regulation rate of ± 2 %. Thus a 
corre ction network is necessary . 


We now refer to the Indicative response curves for the 
values * = 0.1, A = 3.8, f = 0 and M = 1.6. On the curve In 
question we read that the optimized correction network Is a 
network of the type G for which we have a simulation 

value of G - 1.7 and a reduced value for t of 25 (Eq. 3.13). 
If we substitute the true value of G using Eq. (11. 5): 

® " 4.1,7 ^ ■ 0,65 


Since we are using values of M and A which are slightly 
different from those of the simulation, we correct this 
value by using the suggestions given In section 11. A. 1: 


Eq. 


G -:o,es 


1,33 

1,6 


hi 

3.8 


0,65 


We now calculate the value of the 

( 11 . 1 ^ 0 ; 


true time constant using 


7- 25.50. 10~” 



We should therefore Insert the following network: 


v» 

I 

i 


. j 



Gj(p) 


U2.|0*^p • 
2.10“^p 


(12.7) 


We can now refer to the other diagrams concerning this /75 

case and we not that : 

- the overall stability is satisfactory, the reservations 
of a general character expressed on this subject in the 
report on the boost of course remain valid. 

We now see what the minimum power stage gain conditions 

are . 


Since Ig has to vary from 1 A to 0.05 A, we have: 


0,0S 4 J 


( 12 . 8 ) 


Eq. ( 11 . 28 ) gives us, for M = 2 and J = 0.05: 


Hr 5 ^"” 


(12.9) 


The power stage is at its gain minimum for M = 2 and 
A = 22. 


If the correction network was a simple flat gain G = 1.7 
expressed as a simulation value, these new operating conditions 
would not pose any stability problems, proof of which can be 
obtained by looking at the local stability curve. The 
equivalent path band would simply have moved back from. 0.2 P 
to : 

0,2 r - 0,02 F - 400 H« 




But In the Inserted network the zero which is involved 
there is at F/25 » 0.0i| F. There is therefore a risk of 
Instability since the network causes more phase shift at 
0.02 F than at 0.04 F. 

However* given the margin which exists for the indicative 
response diagram of the optimized system, we can move back the 
zero by a factor of three and place it at F/75 = 0.0133 F and 
thus use the network with greater certainty in a vacuum; 


Cj(p) 


0.65 


6.10*^p 


( 12 . 10 ) 


Remark ; If, for example we were working with an instruction 
= 5 V and thus a dividing ridge of ratio 40/5 = 8 on the 
output voltage, it is clear that it would be necessary to 
multiply the gain of the network by 8 in making it change 

to 0.65 X 8 5 . 

12.2 Example 2; Boost Functioning in Continuous 

t, 

Conduction from J=ltoJ=0.2 

Let us take a boost converter having the following 
characteristics: ■ 

Vg - 40 V Xg • I A BMC thus Rg > <0 -^2 

may vary from 30 V to 20 V 

F ■ 20 kHz thus T ■ 50.10”^ 

. . *.4 ( 12 . 11 ) 

We want the converter to operate in ccntlnuous 
conduction from * 1 A up to Ig ^0.2 A, but 
the converter will nevertheless be used from 
Ig - 1 A up to Ig * 0.05 A. There is no 
particular requirement with respect to the 
residual ripple. 

We use a coincidence modulator with A " 5 V. 


This time we will consider Eq. (11.26). The most 
unfavorable case for the transition Into discontinuous 
conduction Is when M - 1.33. Thus this Is In the case when 

we should be at the boundary between the modes. Therefore 
we have : 


M - 1,33 -2m- 2,66 


for J - 0,2 


Hence, for M = 1.33, we have A = 2,66 0.2/1 = 0.53 when J - 1 
We are going to use the following nominal values: 


1*33 

^nc« - °*53 
"sno» - « 


( 12 . 12 ) 


Foi the results relating to the boost we find that It Is 
Inposslbjc to arrive at a satisfactory result with (Ji = 0.1. 
Therefore we assume : 


9 - 0,03 


( 12 . 13 ) 


Eq. (^.3) gives us; 


.40.50.10"* 


,2 0,53 


L ■ 0,7 bB 


( 12 . 114 ) 


Eq. (l4.8) gives us: 


471^. (0,03) 3(2. io4) 2 o,7.10"3 


C - 100 U.F 


( 12 . 15 ) 




I 



The condenser we are using becomes resistive at 10 kHz 
Thus we have: 


n - 2 (12.16) 

The power stage is at its gain maximum for M = 1.33 and 
L " 0«53. 


It is for these values that we are going to calculate the 
correction network. 

Consulting all of the results relative to the boost, we 
then find that for ({> = 0.03 no satisfa ctory solution exis^ 
for a A PS low as 0.53 « 

Two approaches may then be adopted: 

a) We can try to work with a larger A, in other words we 
allow the transition into discontinuous conduction to 
take place for a current clearly greater than 0.2 A. 

This does not seem to us to be a logical solution 
because it amounts to optimizing the converter in the 
continuous mode (which is the most delicate) while it 
will operate most often in discontinuous conduction. 

b) We decrease 4> . For example, we take: H 

(f) = 0.01 (12.17) 

which amounts to using a capacitor: 

C2^^ooO|JLr ( 12 . 18 ) 

In effect, for the results of the boost v:e then find that th£ 
problem has a solution by using a current loop • 
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While observing the development of numeric values when 
i varies from 0.2 to 0.8 then to 2.6. «e see that this 
development is roughly of the type: 




Qj 


fire eoing to interpolate from 
Therefore, as a result, we are going 

the values A = 0.2 and A = 0.8. 

For A = 0.2, we read G = 1.^ and - 1/17. 

For A = 0.8, we read G = 4.6 and G^^ = 1/30. 


Interpolation; 

We have ^ ^ 

4,0 flilY 
thus r^r“^o, 2 y 

476 V O'® / 

17 / 0,8 \* 

likewise 3o"Vo,2y 
. - 0,41 

1/30 0,8 J 


X - 0,858 


0,858 


G ■ 3*2 


, X - - 0,41 




Mote: we may content ourselves with a simplified Interpolation 
in this case. 

NOW let us correct these values to bring them to the 
good valu of M: 


0 - 3.2 




in simulation values 




We calculate the true values: 


G ■ 4.2,66 ^ ■ 1*33 

5*^ 40 - 0,66 (Eq. 11.9) 

When at the pole of the current loop, t = 5 In reduced 
value, thus: 

^ . iJMoli. 4.10"* (En. 11.14) 

27T 

We thus insert the networks: 


G^^p) 

- 1,33 

C»(p) 

0,66 

-5 


U4.10 p 


Finally, we go to the minimum gain conditions of the 
power stage, i .e . M = 2 and J = 0 .5 • 

Using Eq. (11.27) we obtain: 

* I 1 1 , 33 ^ 1 • ,75 

^ »y2,0,53 0^ 2 0,33 

The power stage Is at Its gain minimum when M = 2 and 
when A = 7.5 (for J " 0.05). Therefore, as predicted. It Is 
discontinuous conduction when It will operate correctly for 
the gains in question. 
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As for the elementary units, these reduced parameters will 
characterize each class of PWM shunt converter. 

These parameters are dimensionless and they have the 
greatest possible practical significance. 

For defining these parameters we will make the following 
assumptions; 

- the values of the reduced parameters are defined on the 
basis of the steady state of the system, 

- this steady state is Idealized at the level of the output 
= Vg = cte; Ig = cte. 

These parameters were defined in accordance with the 
ESTEC during the meeting of Jan. 12, 1976. We are going to 
discuss their definition in detail. 

1.1 Reduced Parameter r (GAM) 

It is defined as: r s 2 TT . r .c,.p, 

w * 

F being the operating frequency of the PWM shunt . 

As in the other types of elementary units, this parameter 
characterizes the quality of the output capacitor. If T is 
small, the capacitor is of good quality. If T is large, the 
capcitor Is of poor quality. 

1.2 Reduced Parameter FE (GAE) 

By definition: Pb ^ 2Tr.R.c.r. 


This parameter defines the delay when the diode D of the 
shunt is put into conduction. The greater PE is, the longer the 
delay. 


1.3 Reduced Parameter rS CGAS) ' 


We define; 

Ps = 211'.Rg.C2.F. 

This parameter determines the size of the output capacitor 
C2 of the PWM shunt by comparing the cut-off frequency of the 
output network with the operating frequency of the converter. 

At the same time, this parameter determines the size of the 
residual ripple at the output . 

Parameter MJ, Ratio of Solar Panel Current to Output 

Current 

We set up the following equation; 

A 

Obviously we have : 

0 < M < 1. 

This parameter defines the current operating point of 
the system, l.e. the cyclic ratio of the transistor of the 
power stage. 

1.5 Miscellaneous 

In addition, we have defined the ratio p which exists 
between output resistor Rg and the resistor Rg in parallel on 
the current generator of the solar panel. 



Everything which has already been said in the "Methodology 
Report" in Part I concerning the reduced time constants and the 
gains of the different modulators remains valid. 


Parameters 


in vhat follows we have assumed as known the following: 

- the output voltage Vg, 

- the charging resistance Rg, 

- the operating period of the converter T. 

These quantities Vg, Rg and T relate as well to a 
system as to a real system. We deduce from them the quantiti 
of the circuit to be considered: 

R - 5^-*^ 

Pb.t 


c - 


2vr. 




C2 


IG - 


BC - 


Bg-MJ 

vT 

PS 


TTT nelr, of the Correctlnr. A mpli fier Insuring the Requi red 
± 2 % Regulation 

in the same way as for the other elementary units to 
find out the minimum gain which ensures the 12* regulation, 
suffices to calculate what the maximum variation ,s of the 

control voltage V- of the modulator S 

I vary within the specified limits, then to write the 
s 

following equation: 


For the coincidence modulator we have: 
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- ‘on . ^ 

caax T 




For the magnetic modulators, we have: 


''e-*! " -’■oh- if' 


‘><’*oh <4 


With the definition and the range of variation of MJ 
(cf. section we easily have: 


O.I.T < < 0 . 9 .T 


Thus, finally we have: 


cmin 


- 0.9 A I 

• O.l A J 


coincidence modulator 


Therefore, for the shunt with the coincidence modulator 
we will have: 

0.9 A ^ ^ 

In the case of the simulation, v;e will have: 

SBUOT - 5.6 

For the magnetic modulator we see that we cannot go below 
MJ = 0.333. V -V 

CDOX 1 

V , - 0 

•aln 


V| 

00 - 25 ^-25 because we take Vj - Vj^ 




SBUWr - 25 
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IV. Guidelines for Simulations and Optimizations to be 
Carried Out 
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IV. 1 Variation Range of the Parameters 

In accordance with the ESTEC, we defined on Jan. 12, 1976 
the following values for the reduced pararieters : 


p 

1 

0. 

3. 

7. 

r« 

I 

0.3 

0.6 

1.5 

Ps 

1 

100. 

SCO. 

1000. 

MJ 

t 

0.1 

0.5 

0.9 


From these we deduced the simulations and "optimizations" 
to be carried out which retain a significant value for the 
results . 

IV. 2 Simulations and "Optimizations" to be Carried Out 

- Exhaustive study of the local stability with coincidence 
modulator, 

- simplified study of the local stability with magnetic 
modulator , 

- Investigation of the best correction system for the 
doublets (TE, fS) with the coincidence modulator and 
study of the variations caused by the extreme values of 
MJ and by V , 

- readjustment of the correction system for the magnetic 
modulator , 

- plotting of output Impedances only for optimized systems. 

As for the elementary units studied prevlosuly, we con- 
scructed a flow diagram of the simulations to be carried out. 


All of the doublets 


He . 



All of the doublets 


HJ , n 


ns 



A central value 
of NU 


Local stability with 
coincidence modulator 


r - 0 


Local stability 
with magnetic 
modulator . 


Investigation of 
the best correction 
system with coinci- 
dence modulator 


Two extreme values 
of MJ 

r- 7 


Indicative responses 


^Indicative 

responses 

Output 

impedances 


Adjustment of the 
correction system 


I 

Indicative response with 
magnetic modulator 
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